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Andreev spectroscopy of doped HgTe quantum wells 

M. Guigou and J. Cayssol 

GPMOH, UMR 5798, Umversite Bordeaux I, 

3340,5 Talence, France 

We investigate the Andreev reflection process in high-mobility HgTe/CdTe quantum weUs. We 
find that Andreev conductance probes the dynamics of massive 2+1 Dirac fermions, and that both 
specular Andreev reflection and retroreflection can be realized even in presence of a large mismatch 
between the Fermi wavelengths at the two sides of the normal/superconducting junction. 



I. INTRODUCTION 

Three dimensional topological insulators (TIs) and 
two dimensional Quantum Spin Hall (QSH) states are 
novel electronic phases which are characterized by topo- 
logical invariants rather than by spontaneously broken 
symmetries^"— . Both QSH states and TIs are distin- 
guished from ordinary insulators by the presence of con- 
ducting edge or surface states surrounding an insulat- 
ing bulk'*"^. These protected boundary states were ex- 
perimentally confirmed in HgTe quantum wells^i^, and 
in three dimensional TIs like Bii_3;SbaJ^, Bi2Te3 and 
Bi^Se^iii^. 

All these materials follow a general mechanism 
whereby the strong spin-orbit interaction drives an in- 
version between bands of distinct symmetry^, e.g. oppo- 
site parities^. Simple massive Dirac equations describe 
altogether the conduction/valence and the spin degrees 
of freedom of QSH systems or TIs, the crucial band in- 
version feature being captured by the sign of the mass 
tern*ii"— i^. Therefore doped HgTe quantum wells are also 
model systems to investigate the dynamics of massive 
Dirac fermions in 2+1 dimensions'^"—. 

Combining these unique topological phases with con- 
ventional superconductivity raises fundamental questions 
about topological superconductoraiii^, and may lead to 
the discovery of novel exotic modes such as Majorana 
fermionsi^"— with potential applications for spintron- 
ics and quantum computing^'^ . Previously most stud- 
ies have focused on proximity induced superconductivity 
within ID edge states^i^ or 2D surface statesi^i^"— . 
In doped TIs (resp. QSH systems), the fluctuations 
should be less detrimental to superconductivity than in 
the 2D (or ID) boundary states available in the insu- 
lating regime. Indeed superconductivity was recently re- 
ported in the topological insulator Bi2Se3 doped with 
copper atoms^l while Heusler superconductorsS^ or Thal- 
lium based chalcogenides^ are promising candidates for 
a topological superconductor. 

In this context it is of primary importance to under- 
stand proximity effect between a doped topological in- 
sulator and a s-wave singlet superconductor. The fun- 
damental scattering process is the Andreev reflection 
(AR) whereby an incident electron is converted into a re- 
flected hole at the interface between a normal metal and a 
superconductor—!^. Moreover it was recently predicted 
that relativistic electrons in graphene may experience an 



unusual specular Andreev reflections^ ^. 




(a) 




FIG. 1: (a) Normal/superconductor (NS) contact. (b) 
Schematic band structure where 2|M| is the bulk gap and 
M the mass. The black dots represent the incident and re- 
flected electrons while the open dot stands for the Andreev 
reflected hole (here interband specular reflection is shown). 
Arrows represent group velocities. 



In this paper, we show that AR can be used as a probe 
of the carrier dynamics in HgTe quantum wells which 
is intermediate between linear and quadratic dispersion. 
As low gap semiconductors, HgTe wells also exhibit a 
specular Andreev reflection which is sensitive to band in- 
version. While we use HgTe for calculations, the general 
features are expected to apply for other two-dimensional 
QSH systems"^^ and for three dimensional TIs^'^^'^^. Our 
findings also complete previous studies of the Andreev 
conductance of a NS contact in presence of Rashba spin 
orbital coupling"^^. 

The paper is organized as follows. In Sec II we in- 
troduce the model for the NS contact. The energy and 
angular dependences of the Andreev reflection processes 
are described in Sec. HI. These single channel results 
are used in Sec IV to obtain the multichannel differential 
conductance which is the relevant quantity for transport 
experiments. 



II. MODEL 

In this section, we present the model of the NS junction 
and describe the scattering formalism used to compute 
the Andreev probability and the differential conductance. 



Hamiltonian 



waves with momentum k and energy 



We consider a thin HgTe quantum weU reahzed be- 
tween two identical CdTe barriers. The two-dimensional 
subbands of such quantum wells have been derived from 
the Kane model of HgTe and CdTe using the envelope 
function method^. The QSH state is related to the 
crossing between the electron-like subband El and the 
heavy hole-like subband HI which is controlled by the 
HgTe layer widtb^. When the bulk inversion asymme- 
try is neglected, those relevant bands are degenerated 
with respect to the spin index (±) and have opposite 
parity eigenvalues. Then the low energy dynamics of 
this four-band model is captured by the massive Dirac 
Hamiltonian^ 



_ fh{k) 
^^^)-[ /i*(-k)/' 



(1) 



given in the basis order (|^l-l-) , l-ffl+) , \E1-) , \H1-)). 
The spin up block /i(k) = e(k) -|- (ia(k)(Ta is expressed 
in terms of the standard Pauli matrices CTq (a ~ 1,2,3) 
acting in {\El) , \H1)) space. The Hamiltonian if (k) is a 
Taylor expansion with respect to the in-plane wavevector 
k = {kx, ky) whose coefhcients are constrained by parity 
and time-reversal symmetries. Microscopic theory fur- 
ther yields that da = {Ak:^-Aky,M{k) = M - Bk'^) and 
e(k) = C — Dk^, where k'^ — k^ + ky. The parameters A, 
B, C, D and M depend on the quantum well geometry^. 
In particular, the inversion between \E1) and \H\) is 
controlled by the sign of the mass term M, the QSH 
state being realized in the inverted regime (M < 0). The 
chemical potential C determines the electronic filling of 
the bands which can be electrostatically tuned by the 
action of a distant metallic gate. 

We assume that a singlet s-wave superconductor in- 
duces a pairing potential in some region of a ballistic 
HgTe/CdTe quantum well. Thereby the Hamiltonian of 
Eq. ([T]) must be completed by particle- hole (or charge con- 
jugaison) symmetry leading in principle to an eight by 
eight Bogoliubov-de Gennes Hamiltonian. Since super- 
conductivity only pairs time-reversed states, a spin up 
electron will be coupled with a hole in the spin down 
band. Therefore electrons and holes are described by two 
decoupled (four by four) Bogoliubov-de Gennes-Dirac 
equations like 



h(—ihdr 

A*(r) 



A(r) 
-h{—ihdr 



*(r) = £;*(r). 



(2) 



El- 



* 



Hl- 



and the ab- 



where *(r) = (*£!+, *i/i+, * 
sence of magnetic field is assumed. The excitation energy 
E is measured from the Fermi level. The matrix A is di- 
agonal at the lowest order approximation in momentum. 
In principle the diagonal entries, A^; and A//, of the pair- 
ing matrix are not equal because \E1) and |-ffl) involve 
different combinaisons of atomic orbitals. Nevertheless 
we first assume that A^ — A/z^Aoe'"^ for simplicity. In 
a uniform system, the eigenmodes of Eq.Q are plane 



£;(k) = J(£(k)±d(k))2-fA2, 



(3) 



with d(k) = y^A^k^ + {M -Bk^y. The ± sign refers to 
the conduction/valence subbands. 



B. Scattering problem 

We now consider a straight normal/superconducting 
interface at a; = (Fig[T]). The normal side {x < 0) 
is described by Eq.© with C(r) = Cn and A(r) = 0, 
while the superconducting side may have a distinct elec- 
tronic filling fixed by C(r) = Cs and a uniform pairing 
potential A(r) = AqB^'^. This sharp step model is valid 
when the potentials vary on typical scales smaller than 
the Fermi wavelength of the carriers in HgTe. We fur- 
ther assume a perfect interface which is relevant for high 
quality contacts. Previously a NS contact in the presence 
of Rashba spin-orbit coupling was already considered in 
the context of spintronics^. The main distinction be- 
tween the HgTe/CdTe effective Hamiltonian and similar 
Rashba systems consists in the presence the {m — Bk^)(j3 
term in Eqs. (fTT2l) . 

We solve the scattering problem with an incoming spin 
up quasiparticle in the conduction band. Since the scat- 
tering is elastic and the problem invariant by translation 
along y-axis, the excitation energy E and the transverse 
momentum ky are conserved. Hence all the scattered 
waves are written hereafter as ^(a;)e''^"*'. 

In the normal part, x < 0, electrons and holes are de- 
coupled since A(r) = 0. The electron-like quasiparticles 
are described by four-spinor plane-waves 



*e(a;) 



lx±(ke 



,0,0)e 






(4) 



with x±(k) - (±d(k) -t- Mik),A{kx - iky)}^-2^. The 
upperscript T denotes transpose and the ± sign corre- 
sponds to the conduction or the valence subband. For a 
n-doped HgTe well, the Fermi level lies in the conduction 
subband. The dispersion equation. Cat — £'fc^-|-d(k) = E, 
allows for two possible values of fc^ : 



kUE) = 



2(52-1)2) 



(5) 



± v/72 - 4(i?2 - D'~){M^ - (Cn - E)^)] , 



where 7 = -^2 + 2MB - 2{Cn - E)D and B'^ > D^. 
The positive one, {kf > 0), corresponds to a the prop- 
agative mode in the bulk, while an additional evanescent 
mode (A;| < 0) shows up at interfaces'^"—. The inci- 
dent and reflected electrons have longitidinal momentum 
kex = fcicos^, and — fcicos^ respectively, 6 being the 
incidence angle (FigllJ. The evanescent electron is de- 
scribed by a complex momentum k^x = — *(— ^2 + ^v)^ 



The hole-like quasiparticles are also described by the 
four-spinor plane-waves 



III. SINGLE CHANNEL ANDREEV 
REFLECTION 



^h{x)= (0,0,x±(k/,))e*'='' 



'^"^hx -^ 



(6) 



The longitudinal wavevectors k^x is obtained by solving 
the equation Cat — Dk^ + (i(k) — —E for a hole in the 
conduction band (intraband AR) or Cm — Dk^ — (i(k) = 
— £' for a hole in the valence band (interband AR). 

In the superconducting part {x > 0), the eigenmodes 
are four Bogoliubov quasiparticles which are all evanes- 
cent below the gap {E < Aq) 



*S±(x) = (x±(ks),X±(ks)e*''^)e 



±ip\ {±iksx-K)x 



, (7) 



%.^{x) ^ (x±(k;),X±(k;)e±'^)e(±^'=-±-)-, (8) 



where k = -^/Aq — E"^ /hvp is an inverse coherence 
length. The phase (3 — arccos(i?/Ao) is intrinsically re- 
lated to electron-hole conversion at a normal conductor- 
superconductor interface's. We have introduced the 
wavevectors fcsa; = {klg-kyY^"^ and A:^^ = iiky — klsY^"^, 
where k\g and A:|g are defined by Eq.® under the sub- 
stitution Cn -^ Cs and E — Q. The scattering am- 
plitudes of the four reflected modes and four evanescent 
transmitted modes in the superconductor are determined 
by writing the continuity of the wavefunction and of its 
derivative at x = (Appendix A). 



C. Andreev differential conductance 

When a positive bias V is applied to the normal side 
with respect to the superconductor, the current / is car- 
ried by the injected electrons, the reflected electrons 
(with amplitude ree{E,9)) and the Andreev reflected 
holes (with amplitude rhe{E, 9)). The corresponding dif- 
ferential conductance of the NS interface can be written'^^ 



§ = »<«") 



-i-S 



d9 cos 9 



{1~R{E,9)+Ra{E,9)), 



(9) 



where goieV) = e^ki{eV)W/{Trh), W is the width of the 
HgTe quantum well along y-axis and / = f{E — eV) — 
l/(e^^~'^^-'/"^ -I- 1) is the Fermi distribution of incident 
electrons in the normal lead N at temperature T . The 
probability for an electron of energy E to be reflected as 



an electron is R{E, 9) = 



,{E,9)Y while Ra{E,9) = 



\rhe{E, 



denotes the probability for Andreev 



Ji-Efi) 

reflection as a hole. The average current 

i{E,9) = -2ki{E)cos9{{D + B){d{k)+M{k)f 

+ A^PiD - B) - A^idik) + M{k)}, (10) 

is derived by the standard quantum mechanical proce- 
dure (Appendix B). 



In this section, we investigate the Andreev reflection 
process for a single channel which is labelled by its energy 
E and incidence angle 9. We contrast the behaviors of 
the Andreev reflection probability Ra{E,9) at high and 
low doping levels in the normal part of the NS junction. 
In the later case, the subgap Andreev reflection can be 
either an intraband or interband process whereas it is 
always intraband in heavy doped wells. In both cases, 
varying the ratio B/A allows to turn continuously the 
Andreev reflection probability Ra{E, 9) from the one ob- 
tained in graphene/superconductor junctions to the one 
for standard metal/superconductor junctions. 



A. Heavy doping \Cn\ S> Aq 

We first consider HgTe wells with typical metallic dop- 
ing, i.e. \Cn\ ^ Aq. Then an incident electron (en- 
ergy E) from the conduction band always finds an elec- 
tron of the same band (energy —E) to form a Cooper 
pair thereby realizing an intraband Andreev conversion. 
Therefore the Andreev reflection is a standard retrore- 
fiection because the electron energy E is always much 
smaller than the conduction band Fermi energy \Cn\- 

We find that the probability Ra{E,9 = 0) is strictly 
monotonic (increasing) below the gap, and presents a 
singularity at eV = Aq. The AR probability is con- 
trolled by the Fermi wavelength mismatch (FWM) be- 
tween the normal and superconducting sides, and usually 
decreases when the FWM increases. At given electronic 
fillings (or FWM) , the Andreev probability is suppressed 
when the ratio B/A is increased, namely when going 
from Dirac (purely linear dispersion) to Schrodinger 
(purely quadratic dispersion) dynamics (Fig. [2]). This 
crossover is characterized by the dimensionless param- 
eter A = ki{0)B/A= (n/no)^/^ where n = fc?(0)/(27r) 
is the two-dimensional carrier density in the conduction 
subband while no = A^/{2nB'^) ~ S.IO^^ ^^-^-2 for A = 4 
eV • A and B = —70 eV • A^. When electrons are injected 
through the junction with a finite incidence angle, the 
Andreev reflection probability decreases but has qualita- 
tively the same energy dependance than the 9 = mode 
(Fig. [5]). The general crossover between purely linear 
or quadratic dispersion still pertains for any incidence 
angle. 



B. Low doping \Cn\ < Aq 

We now turn to the case of low doping in the nor- 
mal side, namely \Cn\ < Aq. At low voltage bias 
(0 < eV < \Cn\ - \M\), the reflected hole belongs to 
conduction band thereby realizing the standard intra- 
band AR which is a retroreflection. At intermediate 
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FIG. 2: (Color Online). Single channel Andreev reflection 
probability through the NS junction for 5 = (top panel) 
and 6^ = 7r/10 (bottom panel). The electronic fillings are set 
by — Cjv = 30 meV and —Cs = 1 eV, while the induced gap is 
Ao = 1 meV. All curves are plotted for A = 4 eV-A, and in the 
inverted regime A/ — — O.IAq. The Andreev reflection proba- 
bility is progressively suppressed when the quadratic disper- 
sion is increased: -B = 0.7; 7; 70; 700 and 7000 eV ■ A^ 
from top to bottom. In contrast to B, the parameter D has 
very little impact on Andreev reflection probability provided 
\B\ > \D\. Hence we have set D = in order to allow for a 
broader window of variation for the parameter B. 
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FIG. 3: (Color Online). Single channel Andreev reflection 
probability in the inverted regime (M = — O.lAo) at low 
doping (—Cat = 0.5Ao) and for three different incidence an- 
gles: — Q (blue solid line), 9 = 7r/30 (red dashed line) 
and 6 — tt/IO (green dotted line). The other parameters are 
taken from tyisical HgTe wells: A = 4 eV- A, B = -70 eV- A^ 
D = —50 eV-A^. The superconductor is heavily doped (Fermi 
energy —Cs = 1 eV). 



reflection is totally suppressed {Ra{E,9) = 0) broad- 
ens while increasing incidence angle (Fig. [2|). This phe- 
nomenon has the following explanation. Owing to trans- 
lational invariance along the NS interface, the transverse 
momentum ky is conserved and one has to solve an effec- 
tive one-dimensional scattering problem for each value of 
ky. It turns out that a finite ky acts as a supplementary 
mass/gap in this ID scattering problem. For instance, 
when M = the 2D system is gapless but the ID prob- 
lem (described by the dispersion relation E{kx,ky) for 
defined ky) acquires a gap proportionnal to the trans- 
verse momentum, as it is well known in graphene^. 



IV. DIFFERENTIAL CONDUCTANCE 



vohages (ICatI - \M\ < eV < \Cn\ + \M\), the AR 
is completely suppressed in the bulk because the hole 
should be in the semiconducting gap. At higher volt- 
ages {eV > I Cjv I + l-^-^l), the reflected hole corresponds 
to the removal of a valence band electron. Then the 
AR is an interband process and a specular reflection. In 
brief, Andreev reflection allows for a spectroscopy of the 
HgTe/CdTe well 2D spectrum since for instance the en- 
ergy window where the AR is totally suppressed is 2 \M\, 
i.e. twice the semiconducting gap, for the mode 6 = 0. 

Moreover at the bottom of the conduction band, a sin- 
gularity appears in the AR probability at _E = | Cn | — | ^^ | 
and normal incidence (blue solid curve in Fig. [3]). Nev- 
ertheless, this abrupt behavior becomes smoother when 
electrons are injected with an non zero angle and the AR 
probability decreases. 

Furthermore the energy window where the Andreev 



In this section, we focus on the differential dl/dV con- 
ductance of the NS contact which is obtained by angular 
integration over all the transverse channels. The dl/dV 
characteristics are qualitatively distinct in the low and 
in the heavily doped regimes respectively, and they de- 
pend quantitatively upon the ratio B/A and the Fermi 
wavelength mismatch (FWM). 



A. Heavy doping |Cjv| ^ Ao 

The voltage dependance of the differential conductance 
dl/dV is inherited from the energy dependence of the 
AR probabflity through Eq.dH]). Hence, in the case of 
heavy doping in the N region, the dl/dV characteristics 
are also strictly increasing below the gap, and present 
a singularity at eV = Aq (compare Fig. H] and Fig. [5]) . 
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FIG. 4: (Color online). Differential conductance of NS con- 
tact normalised by go{eV) = e^ki{eV)W/{'Kh). Both side of 
the junction are heavily doped —Cn = 30 meV and —Cs = 1 
eV. All curves are plotted for A — 4 eV • A, and in the in- 
verted regime M — — O.lAo. The Andreev conductance is 
plotted for differents values of B : ~B = 0.7; 7; 70; 700 and 
7000 eV • A^ from top to bottom. As in Figs. [2l the red 
dotted curve is identical to the graphene case up to a factor 
2 due to the absence of valley degeneracy in HgTe and the 
blue solid curve corresponds to typical parameters for current 
HgTe/CdTe. We have set D = 0. 



The differential conductance of a normal potential step is 
asymptotically recovered far above the gap {eV ^ Aq). 

We now discuss the crossover between Dirac (purely 
linear dispersion) and Schrodinger (purely quadratic dis- 
persion) dynamics (Fig. \^. For this purpose we use 
again the dimensionless parameter A — ki{0)B/A de- 
fined in section III. A. The increase of the ratio B/A 
leads to the decay of the Andreev conductance. This re- 
sult is in accordance with the fact that FWM strongly 
suppresses AR in standard metalg'^^, whereas AR is 
very robust against FWM at a graphene/superconductor 
interface^i^. These contrasted behaviors can be used to 
probe the relative strength of linear and quadratic dis- 
persion in HgTe wells (Fig. |4]). 

Finally thin films of 3D TIs are ideal systems to probe 
the complete crossover regime since their effective pa- 
rameters A and B can be tuned by varying the film 
thickness^, whereas only Af can be significantly varied 
in HgTe quantum wells. 



cess (0 < eV < \Cn\ ~ \M\), forbidden in the bulk 
{{\Cn\ - \M\ < eV < \Cn\ + \M\)) and an interband 
process {eV > \Cn\ -f |M |) which is reminiscent of spec- 
ular AR in grapheno22ii2^. Unfortunately the observation 
of specular AR in graphene was hindered so far by disor- 
der effects occuring at low doping. We expect that spec- 
ular AR should be more easily observed in HgTe quan- 
tum wells since actual samples are characterized by larger 
elastic mean free paths than in graphene on a substrate. 
Unlike the Andreev reflection probability, the behav- 
ior of the differential conductance in the inverted regime 
does not emphasize any peak at energy corresponding to 
the bottom of the conduction band. At such energy, only 
one channel participates to the Andreev conductance so 
that its contribution is tiny compared to the number of 
channels that are considered in the the ballistic conduc- 
tance go{E). 




B. 



Low doping |Cjv| < Ao 



FIG. 5: (Color online). Differential conductance in the non 
inverted (Af — O.lAo, dashed red line) and inverted (M = 
— O.lAo, solid blue line) regime with —Cn = 0.5Ao. The 
parameters A, B and D are the same than in Fig. [S] 

Moreover there is a sizeable difference between the nor- 
mal {A4/B < 0) and the inverted regimes [M / B > 0) 
(Fig. [5]). Below the gap, the interband specular AR is 
enhanced in the non inverted regime with respect to the 
inverted regime, whereas Andreev conduction is stronger 
in the inverted regime when < eV < \Cn\ — \M\ (in- 
traband AR) or above the gap. It should be empha- 
sized that, even at the band crossing M — 0, the AR in 
HgTe/CdTe quantum wells differ from graphene owing 
to the presence of quadratic terms Bk^ and Dk^ in the 
Hamiltonian Eq.([5]). 



In the case of low doping in the normal side {\Cn\ < 
Ao), the subgap differential conductance reveals both the 
semiconducting gap and the sign of the mass M (Fig. 
[5]). Foremost the behavior of the differential conduc- 
tance shows the presence of the reffected hole in the 
conduction/valence band through the electron-hole con- 
version at the interface. As a function of the voltage 
bias, the AR process is successively an intraband pro- 



C. Experimental realization 

The predicted differential conductance can be checked 
experimentally by realizing transparent enough contacts 
on existing mercury telluride heterostructures. In partic- 
ular at low doping, the specular (interband) Andreev re- 
flection could be observed in HgTe/CdTe quantum wells 
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FIG. 6: (Color online). Influence of the Fermi energy ~Cs in 
the superconductor while — Cjv = 0.5Ao, M — — O.IAq. The 
values of the parameters A, B and D are the same than in 
Fig. El 



with a superconducting contact. Nevertheless mass fluc- 
tuations, related to spatial variations of the HgTe layer 
width, may hinder the observation of specular AR in such 
low doped HgTe/CdTe wells (Fig. [5]). 

Another difficulty is related to the fact that the in- 
fluence of doping in the superconducting region might be 
difficult to separate from the one of other parameters, like 
the mass M or the ratio B/A, using Andreev measure- 
ments only (Fig. [S]). Moreover \Cs\ is neither controlled 
nor even known in contrast to \C'n\ which can be tuned by 
gating. We anticipate that future experiments will solve 
this issue by combining transport measurements through 
NN and NS junctions. 

Finally, one should also expect in the inverted regime 
a contribution from the helical edge states. Since 
there is perfect Andreev reflection at a QSH insula- 
tor/superconductor junction, this contribution should be 
4e^//i (2e^//i per edge)SS. In the metallic regime, the 
large number of channels makes this edge conductance 
negligible with respect to the bulk conductance. Nev- 
ertheless near the band edges, i.e. at very low doping, 
the edge contribution becomes comparable with the bulk 



V. CONCLUSION 

In conclusion, Andreev conduction is a probe of the 
carrier dynamics in doped HgTe/CdTe quantum wells. 
Indeed the underlaying AR mechanism is extremelly 
sensitive to the balance between Schrodinger (purely 
quadratic dispersion) and Dirac (purely linear disper- 
sion) dynamics. The unavoidable FWM between the 
normal and superconducting sides is not so harmful in 
HgTe/CdTe wells than in usual metals (but also more 
detrimental than in graphene). Furthermore we ex- 
pect that our analysis of intraband AR pertains for 



strong TIs in the doped regime, like Bi2Se3 or Bi2Te3, 
since they are described by similar effective Hamiltonians 
than HgTe/CdTe quantum wellsiii. Nevertheless those 
large gap 3D TIs are less favorable than the small gap 
HgTe/CdTe wells for the observation of specular AR. 
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VI. APPENDIX 
A. Boundary conditions 

We provide a derivation of the boundary conditions 
we used in this paper to determine the scattering ampli- 
tudes. 

The Bogoliubov-de Gennes matrix equation Eq.([5]) 
consists in four scalar equations. For definiteness let us 
consider the first scalar equation: 

{C + A^Ei+ix, y) + {D + B){dl + dl)^Ei+{x, y) 



-iAd:,-^ Hi+{x,v) + Ady-^ Hi+{x,v) 
+e(a:)Aoe**n^_(a;,y) = E-^ei+{x,v), 



(11) 



where the Heaviside function Q{x) indicates the presence 
of superconductivity only on the right side of the junc- 
tion and where we made the substitution kx^y — >■ —idx,y 
Integrating Eq. ljlip over the small region [— e; e] around 
the interface, we obtain 



{C + M) 



dx^ 



El- 



Ax, y) 



(12) 



{D + B)j dx(dl+dl)-^Ei+{x,v) 

dxdx'iHi+{x,y)+ Ady dx-^Hi+{x,y) 



-iA 



+Aoe*'^ / dx^*E^_^{x,y) ^E dx^Ei+{x,y), 

Taking into account that the spinor ^(a;, y) = 
(*Ei+(a;,y),*ffi+(x,2/),*|5i_(a;,y),*^i_(a;,y)) is 

bounded in the region [— e;e], the limit e — >■ yields the 
relation 

lim{D + B) [a,*£i+(6, y) - a,*£i+(-e, y)] 

-limzA['J'Hi+(e,y)-*m+(-e,2/)] =0. (13) 

The second term of the left-hand side is equal to zero, 
according to the continuity of the wavefunctions (other- 
wise the first order derivative are not defined). Hence 
dx^ Ei+ix^y) is also continuous at a; = 0. 



Following the same procedure for the three others dif- 
ferential equations, we obtain the complete set of bound- 
ary conditions 






B. Average currents 



<^fr — 



dH 

dkx 



I 


-2D+kx 


A 







A 


-2L»_fc^ 










2D+kx 


-A 


\ 





-A 


2D^k 



(14) 



We give here the derivation of the average current. The 
current operator in the a;-direction is defined bjJ^ 



where H is the Bogoliubov-de Gennes Hamiltonian ap- 
pearing in Eq.(I2|) and with D± = D ± B. The average 
current for spin-up electrons writes 



Je{E) = *:(x)J.*e(a 



(16) 



where *e(a;) is defined in Eq.(|4l). Substituing Eas. (|4ll5p 
into Eq. P^ leads to the expression Eq. pUj) of the aver- 
age current. 

The average current for a reflected hole can be obtained 
by following the same procedure under the substitutions 
\l/e(x) — >■ ^h{x) and kex — >■ khx- Then the average cur- 
rent of refiected holes is 

jhiE,e) = T2ki{~E)cose{D+{d{k) + M{k))^ 

+ A^k^D^ - A2(d(k) + M(k)}, (17) 



(15) where the =F sign refers to the conduction/valence band. 
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